DTiC  file  COW 


AGARD-AR-267 


DISTRIBUTION  AND  AVAILABILITY 
ON  BACK  COVIR 


'Mpprarid  tar  pRUta  _ 

xmmrrmmr  tuMM 


«*** 


90  06  27  133 


AGARD-AR-267 


NORTH  ATLANTIC  TREATY  ORGANIZATION 
ADVISORY  GROUP  FOR  AEROSPACE  RESEARCH  AND  DEVELOPMENT 
(ORGANISATION  DU  TRAITE  DE  L'ATLANTIQUE  NORD) 


AGARD  Advisory  Report  No.267 

TECHNICAL  EVALUATION  REPORT 
on  the 

FLUID  DYNAMICS  PANEL  SYMPOSIUM 
on 


Computational  Methods  for  Aerodynamic  Design 
(Inverse)  and  Optimization 


(Les  methodes  de  calcul  pour  la  conception  aerodynamique 
(methodes  inverses)  et  l’optimisation) 
by 

Preston  A.Hennc 

Manager  —  MD-90  Technical  Integration 
Douglas  Aircraft  Company 
McDonnell  Douglas  Corporation 
3855  Lakewood  Blvd 
Long  Beach,  Ca  90846 
United  States 

Edited  by 

Professor  Ir  J.W.Slooff 
National  Aerospace  Laboratory,  NLR 
1059  Amsterdam 
Netherlands 


This  Advisory  Report  was  produced  at  the  request  of  the  Fluid  Dynamics  Panel  of  AGARD, 


The  Mission  of  AGARD 


According  to  its  Charter,  the  mission  of  AGARD  is  to  bring  together  the  leading  personalities  of  the  NATO  nations  in  the 
fields  of  science  and  technology  relating  to  aerospace  for  the  following  purposes: 

—  Recommending  effective  ways  for  the  member  nations  to  use  their  research  and  development  capabilities  for  the 
common  benefit  of  the  NATO  community; 

—  Providing  scientific  and  technical  advice  and  assistance  to  the  Military  Committee  in  the  field  of  aerospace  research 
and  development  (with  particular  regard  to  its  military  application); 

—  Continuously  stimulating  advances  in  the  aerospace  sciences  relevant  to  strengthening  the  common  defence  posture; 

—  Improving  the  co-operation  among  member  nations  in  aerospace  research  and  development; 

—  Exchange  of  scientific  and  technical  information; 

—  Providing  assistance  to  member  nations  for  the  purpose  of  increasing  their  scientific  and  technical  potential; 

—  Rendering  scientific  and  technical  assistance,  as  requested,  to  other  NATO  bodies  and  to  member  nations  in 
connection  with  research  and  development  problems  in  the  aerospace  field. 

The  highest  authority  within  AGARD  is  the  National  Delegates  Board  consisting  of  officially  appointed  senior 
representatives  from  each  member  nation.  The  mission  of  AGARD  is  carried  out  through  the  Panels  which  are  composed  of 
experts  appointed  by  the  National  Delegates,  the  Consultant  and  Exchange  Programme  and  the  Aerospace  Applications 
Studies  Programme.  The  results  of  AGARD  work  are  reported  to  the  member  nations  and  the  NATO  Authorities  through  the 
AGARD  series  of  publications  of  which  this  is  one. 

Participation  in  AGARD  activities  is  by  invitation  only  and  is  normally  limited  to  citizens  of  the  NATO  nations. 


The  content  of  this  publication  has  been  reproduced 
directly  from  material  supplied  by  AGARD  or  the  authors. 


Published  April  1990 

Copyright  C  AGARD  1990 
All  Rights  Reserved 

ISBN  92-835-0557-3 


Primed  by  Specialised  Priming  Services  Limited 
40  Chigweli  Lane,  Loughton,  Essex  IG/0  JTZ 


Foreword 


The  Fluid  Dynamics  Panel  of  AGARD  organised  a  Specialists'  Meeting  on  the  subject  of  "Computational  Methods  for 
Aerodynamic  Design  (Inverse)  and  Optimization".  The  Specialists’  Meeting  was  motivated  by  the  observation  that  “design 
type"  of  CFD  (Computational  Fluid  Dynamics)  methods  appear  to  receive  relatively  little  attention  as  compared  to  "analysis 
type"  of  methods;  this  in  spite  of  the  fact  that  the  “design  type"  of  methods  offer  unique  possibilities  for  which  there  is  no 
equivalent  in  experimental  aerodynamics. 

The  Program  Committee  for  the  meeting  is  grateful  that  Mr  Preston  Hcnne  of  McDonnell  Douglas  Corporation,  Douglas 
Aircraft  Company  accepted  the  invitation  for  acting  as  the  Technical  Evaluator  of  the  meeting,  in  particular  because  he 
combines  a  detailed  knowledge  of  "design  type”  CFD  methods  with  a  vast  experience  in  aerodynamic  design  in  an  industry 
environment.  The  present  report  contains  his  observations,  remarks  and  comments  on  the  meeting.  The  23  papers  presented  at 
the  meeting  have  been  collected  in  AGARD  Conference  Proceedings  CP  463. 


Le  Panel  AGARD  de  la  Dynamique  des  Fluides  a  organise  une  reunion  de  Specialistes  sur  “Les  methodes  de  calcul  pour 
la  conception  aerodynamique  (methodes  inverses)  et  I’optimisation”.  En  effet.  le  Panel  a  constate  que  les  methodes  CFD 
(I'aerodynamique  numerique)  pour  la  conception  suscitent  relativement  peu  d  interet  par  rapport  aux  methodes  pour  I'analysc. 
et  ceci  malgre  le  fait  que  les  methodes  pour  la  conception  offrent  des  possibilites  uniques,  pour  lesqucllcs  il  n’existe  aucun 
equivalent  dans  le  domaine  de  I’aerodynamique  experimentale. 

Le  comite  du  programme  dc  cette  reunion  tient  a  remercier  M.  Preston  Henne  de  la  McDonnell  Douglas  Corporation. 
Douglas  Aircraft  Company,  pour  avoir  bien  voulu  acccpte  d’exercer  les  functions  d’Expert  technique  pour  cette  reunion.  Le 
comite  se  felicite  sur  son  choix  dexpert.  puisquil  s’agit  de  quelqu’un  qui  sait  allier  des  connaissances  approfondies  des 
methodes  CFD  du  type  “conception"  a  une  vaste  experience  dans  le  domaine  de  la  conception  aerodynamique.  dans  un 
contexte  industriel. 

Ce  rapport  contient  ses  reflexions,  remarques  et  commentaires  sur  la  reunion.  Les  23  communications  presentees  lorsde 
la  reunion  ont  ete  rassemblees  sous  la  forme  du  Compte-rendu  de  Conference  AGARD  CP  463. 


J.W.SIooff 
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AGARD  FLUID  DYNAMICS  PANEL  -  64th  MEETING 

SPECIALISTS'  MEETING  ON 

COMPUTATIONAL  METHODS  FOP  AERODYNAMIC  DESIGN 
( INVERSE)  AND  OPTIMIZATION 

LOEN,  NORWAY 
22-23  MAY  1989 


by 

Preston  A.  Henne 

Manager  -  MD-90  Technical  Integration 
Douglas  Aircraft  Company 
McDonnell  Douglas  Corporation 


3855  Lakewood  Blvd. 
Long  Beach,  Ca.  90846  USA 


1.0  SUMMARY 

The  papers  presented  at  the  AGARD  Specialists'  Meeting  on  Computational  Methods  for  Aerodynamic  Design 
(Inverse)  and  Optimization  are  reviewed.  Strengths  and  weaknesses  are  Identified  for  many  of  the 
contributions  as  each  is  reviewed.  The  reviewer  closes  with  some  general  comments. 


2.0  INTRODUCTION 

The  64th  Meeting  of  the  AGARD  fluid  Dynamics  Panel  Included  a  Specialists'  Meeting  on  22-23  May, 
1989.  The  Specialists'  Meeting  was  entitled  "Computational  Methods  for  Aerodynamic  Design  (Inverse)  and 
Optimization."  The  Program  committee.  Appendix  A,  Identified  the  theme  for  the  meeting  as  the  following: 

"Computational  Fluid  Oynamlcs  (CFD)  play  an  Increasingly  Important  role  In  aerodynamic  design. 
From  the  design  applications  point  of  view  two  catergorles  of  CFD-based  design  methodology  may  be 
distinguished. 

The  first  utilizes  analysis  type  methods  In  an  heuristic/empirical  cut-and-try  type  of  process. 
In  this  kind  of  process  the  role  of  CFO  Is  to  predict  the  aerodynamic  characteristics  of  a 
configuration  (or  part  thereof)  of  given  geometry  provided  by  the  designer.  The  second  category  of 
CFD-based  methods  addresses  the  problem  of  design  for  aerodynamic  characteristics  In  a  more  direct 
sense.  Examples  are  (the  classical)  Inverse  methods  which  provide  the  detailed  geometry  required  to 
generate  a  given  pressure  distribution  and  methods  utilizing  numerical  optimization  techniques  to 
obtain  the  geometry  that  minimizes,  subject  to  constraints,  a  given  aerodynamic  objective  function 
such  as  drag,  load  distribution,  etc. 

The  purpose  of  the  Specialists'  Meeting  Is  to  stimulate  communication  on  recent  developments  and 
current  research  on  the  second  category  of  methods  (f.e.  Inverse  methods  and  methods  utilizing 
numerical  optimization  techniques).* 

The  Specialists'  Meeting  was  organized  Into  the  following  four  sessions: 

Session  I  -  Invited  and  Survey  Papers 

Session  II  -  Inverse  Methods  /  Airfoils  and  Mings 

Session  III-  Inverse  Methods  /  Turbomachinery,  Intakes,  Ducts 

Session  IV  -  Numerical  Optimization 

The  full  listing  of  authors  and  papers  Is  Included  In  Appendix  B.  The  paper  number,  shown  In  Appendix  8 
and  utilized  In  the  discussion  below.  Is  the  number  assigned  by  the  program  committee  and  does  not 
represent  the  order  of  the  program  presentations. 


3.0  SESSION  I  -  INVITED  AND  SURVEY  PAPERS 

In  thl*  introductory  session  four  presentations  were  made.  The  papers  of  Sobleczky,  Boccl,  and  Koster 
et.  al .  were  available  at  the  meeting.  The  fourth  presentation,  that  of  Jameson,  was  an  oral  briefing 
only.  The  focus  of  this  Introductory  session  was  largely  on  airfoils  and  wings.  The  technical  range  and 


contrast  of  the  papers  In  this  session  was  quite  large.  The  presentations  Included  reviews  of  remarkably 
elegant  and  new  design  methods  as  well  as  applications  of  traditional  methods  with  reportedly 
disappointing  results.  The  variation  In  success  reported  In  these  Introductory  papers  was  larger  than 
expected. 

PAPER  1.  SOBIECZKY  provided  the  extended  Introductory  presentation  and  focused  on  flowfleld 
characteristics,  particularly  as  Inferred  from  hodograph  plane  analyses.  The  printed  version  of  the 
Sobleczky  paper  highlights  the  different  design  method  approaches.  These  approaches  are  the  hodograph 
method,  the  Inverse  method,  and  parametric  or  numerical  optimization.  The  written  version  of  the  paper 
contains  a  brief  summary  of  recent  work  In  each  category  of  design  method  approaches. 

In  the  Inverse  method  area  In  particular,  the  Euler  equation  airfoil  design  method  of  Drelal  was 
highlighted.  The  Orela  method  Is  a  novel  approach  to  airfoil  design  that  Is  a  natural  extension  of  the 
streamline  oriented  numerical  scheme  originally  developed  for  an  analysis  method.  The  work  of 
TakanashK  is  also  mentioned.  Takanashl  utilized  an  approach  based  on  a  residual-correction  scheme 
coupled  to  transonic  Integral  equation.  The  advantage  of  this  approach  Is  that  the  geometry  correction 
scheme  can  be  maintained  separate  and  distinct  from  the  flow  solver.  The  flow  solver  Is  treated  as  a 
■black  box".  Hence,  the  scheme  Is  applicable  In  principle  to  many  existing  analysis  methods. 

Numerical  optimization  works  of  Consentlno  and  Holst3  and  of  Gregg  and  Mlsegades*  are  mentioned. 
Consentlno  and  Holst  utilized  a  gradient  search  strategy  In  combination  with  a  parameterized  definition  of 
a  portion  of  a  wing  surface.  The  objective  function  for  the  numerical  optimization  was  1 1 ft-to-drag. 
Gregg  and  Mlsegades  performed  similar  wing  optimization  but  utilizing  an  evolution  theory  numerical  scheme 
In  place  of  a  gradient  based  method.  The  advantage  reported  for  the  evolution  theory  scheme  Is  the 
ability  to  handle  a  large  number  of  design  variables,  coupled  design  variables,  complex  constraints,  and 
step  functions. 

Sobleczky  closes  his  paper  with  a  forecast  dealing  with  Artificial  Intelligence  and  Expert  Systems. 
Sobleczky  states  that  progress  should  be  expected  In  the  development  of  Aerodynamic  Expert  Oeslgn  Systems 
In  the  near  future. 

PAPER  3.  BOCCI  presented  a  review  of  airfoil  design  techniques  and  design  exercises  accomplished  at 
the  Aircraft  Research  Association  Limited  (ARA)  over  a  number  of  years.  Oeslgn  studies  of  a  combat  wing 
airfoil,  a  transport  wing  airfoil,  laminar  flow  airfoils,  and  propeller  blade  airfoils  are  reviewed  In  a 
summary  fashion. 

This  paper,  different  from  any  of  the  other  presentations  made  at  the  conference,  seemed  to  be  quite 
negative  on  the  use  the  recently  developed  design  methods.  Examples  are  shown  which  are  reported  to  be 
failures  of  the  design  methodology,  particularly  Inverse  supercritical  methods.  However,  In  light  of  a 
growing  volume  of  work  to  the  contrary  and  reported  by  user  agencies  (see  References  5-10  as  a  quick 
sample),  this  position  seems  difficult  to  understand.  Boccl's  paper  seems  out  of  line  with  current 
capabilities  and  the  state-of-the-art  In  design  methodology.  Review  of  the  other  papers  at  this 
conference  also  support  such  a  view. 

PAPER  7.  KOSTER  et.al.  provided  only  hard  copies  of  vlewfoll  art  at  the  time  of  the  conference. 
Hence,  It  Is  somewhat  more  difficult  to  critique  such  an  effort.  In  the  presentation  reference  was  made  to 
at  least  five  different  design  methods  used  at  the  DRl-Instltute  for  Design  Aerodynamics,  FRG.  This 
reference  Indicates  a  strong  commitment  to  design  methods  at  DRL .  Examples  were  presented  for  an  airfoil 
application  and  a  nacelle  cowl  development.  Koster  makes  the  point  that  four  essential  parts  are  needed 
for  the  design  process:  (1)  detailed  and  accurate  description  of  the  design  requirements,  (2)  design 
methods  for  providing  the  basic  shape  and  for  carrying  out  small  changes,  (3)  analysis  methods  to  confirm 
the  design  and  estimate  .he  off-design  behavior,  and  (4)  experience  In  the  use  and  combination  of  results 
from  design  and  analysis  methods  In  order  to  perform  a  successful  design.  The  fourth  Item  cannot  be 
overemphasized.  The  finest  methods  available  today  still  demand  application  experience  If  they  are  to  be 
used  with  confidence. 

PAPER  22.  JAMESON  presented  an  oral  review  of  recent  work  that  he  has  accomplished  In  the  area  of 
aerodynamic  design.  Specifically,  Jameson  has  developed  a  design  method  approach  by  applying  control 
theory  to  the  problem.  He  treats  the  design  problem  as  a  control  problem  In  which  the  control  Is  the 
shape  of  the  aerodynamic  surface.  By  using  control  theory  a  target  pressure  distribution  can  be  sought 
while  additional  quantities  such  as  drag  can  be  minimized  simultaneously. 

The  elegant  mathematical  formulation  presented  by  Jameson  Is  clearly  the  product  of  some  creative 
thinking  focused  on  the  the  aerodynamic  design  problem.  Jameson  Indicated  that  such  a  scheme  has  been 
mathematically  derived  for  two-dimensional  potential  flow,  two- d Imes Iona 1  Euler  flow,  and 
three-dimensional  Euler  flow.  Computational  Implementation  of  the  scheme  has  been  accomplished  for  the 
two-dimensional  flows.  Examples  were  shown  that  Illustrate  a  significant  potential  for  the  method.  A 
drawback  of  the  scheme  as  It  currently  stands  Is  the  computational  time  Involved.  Each  Iteration  of  the 
scheme  requires  about  the  time  of  two  flow  solutions.  A  second  differential  equation,  the  adjoint 
equation  with  Its  related  boundary  conditions,  must  be  solved  along  with  a  standard  potential  equation  or 
Euler  equation  solution  In  each  geometry  Iteration.  If  5  to  10  geometry  iterations  are  required  for 
convergence,  then  the  process  Is  roughly  equivalent  to  10  to  20  analysis  solutions.  Nevertheless,  this 
approach  represents  a  fresh  look  at  the  design  problem  and  numerous  researchers  are  likely  to  adopt  It. 


4.0  SESSION  II  -  INVERSE  METHODS  /  AIRFOILS  AND  WINGS 


PAPER  2.  VOLPE  discussed  recent  work  assisted  with  transonic  airfoil  design.  The  theme  of  his 
presentation  was  based  on  highlighting  various  approaches  to  satisfying  Llghthl 1 1 ' s  three  constraints. 
The  three  constraints  are  associated  with  freestream  consistency  and  orthogonal  trailing  edge  closure  of 


3 


the  designed  airfoil .  Volpe  has  shown  that  these  constraints  can  be  Introduced  Into  the  design  problem  In 
an  number  of  ways  and  has  used  the  Olrlchlet  boundary  condition  (transpiration  scheme)  In  a  2-0  potential 
solution  to  demonstrate  the  approach.  A  number  of  sample  applications,  both  subsonic  and  transonic,  of 
the  design  scheme  were  presented.  Oetalls  were  presented  that  substantiate  the  method  as  being  a  well 
developed  and  matured  capability. 


PAPER  ♦.  MALONE  presented  results  for  an  application  of  a  residual-correction  approach  coupled  to  a 
2-0  Navler-Strokes  code.  As  discussed  In  the  written  version  of  the  paper,  this  approach  allows  the 
geometry  correction  scheme  to  drive  the  analysis  code  as  a  “black  box".  The  title  of  the  paper  Is  "An 
Efficient  Airfoil  Design  method  Using  the  Navler-Stokes  Equations".  Clearly  this  work  Is  one  of  the 
earliest  attempts  at  using  the  N-S  equations  In  an  Inverse  method  application.  Reported  solutions 
utilized  several  thousand  flowfleld  Iterations  for  convergence.  This  magnitude  of  resource  use  Indicates 
a  different  definition  of  efficiency  when  compared  to  that  associated  with  potential  solution  schemes. 

Three  example  applications  are  shown  In  the  paper.  The  first  two  demonstrate  solution  convergence 
while  the  third  represents  a  more  typical  design  application.  Examination  of  the  results  for  the  third 
application  reveals  that  a  target  pressure  distribution,  defined  parametrically  a  priori,  was  closely 
achieved  In  the  design  process.  However,  the  resultli..  airfoil  shape  would  appear  to  have  an  unusual 
shape  near  the  leading  edge  and  a  surface  wiggle  at  the  shock  position  on  the  upper  surface. 

It  was  the  reviewer's  hope  that  the  author  had  actually  attempted  a  design  at  massively  separated  flow 
conditions,  thus  Justifying  the  use  of  the  N-S  set  of  equations.  Unfortunately,  this  was  not  the  case  and 
only  an  attached  flow  case  was  presented  at  this  time.  Future  efforts  would  be  expected  for  separated 
flow  conditions  In  order  to  Justify  N-S  application. 

PAPER  5.  FORNASIER  described  a  progress  report  on  Inverse  capability  using  a  higher-order  panel 
method  for  arbitrary  aircraft  configurations  at  subsonic  and  supersonic  speeds.  The  panel  method  Is  based 
on  mixed  boundary  conditions  of  Olrlchlet  and  Neumann  type.  Inverse  capability  Is  Introduced  by  defining 
additional  singularity  distributions  In  an  Iterative  fashion  to  drive  an  Initial  configuration  to  a  new 
configuration  with  a  specified  pressure  distribution.  Examples  are  presented  In  the  written  paper  for  a 
2-0  airfoil  case  and  for  a  3-0  body  case.  In  both  cases  convergence  toward  the  correct  geometry  for  a 
corresponding  target  pressure  distribution  Is  demonstrated.  Detailed  review  of  the  results  still  pointed 
out  some  modest  discrepancies,  according  to  the  author,  and  future  efforts  are  directed  at  further 
generalizing  the  geometry  handling  procedures. 

PAPER  6.  BRANOSHA  presented  a  scheme  for  transonic  wing  design  based  on  a  residual-correction 
approach.  In  this  scheme  an  Initial  3-0  wing  geometry  Is  used  to  create  an  Initial  pressure  distribution 
using  the  transonic  analysis  code  XFL022.  A  target  pressure  distribution  Is  specified  and  compared  to  the 
Initial  pressure  distribution.  This  comparison  yields  a  pressure  difference  referred  to  as  a  defect 
pressure.  This  Incremental  Cp  Is  transformed  Into  an  equivalent  subsonic  Incremental  velocity  which  Is 
Imposed  on  a  3-D  subsonic  Inverse  panel  method  based  on  mean  plane  singularity  distributions.  Further 
constraints  can  be  Imposed  on  camber  and  thickness  and  a  least  squares  solution  Is  obtained  to  minimize 
the  pressure  defect  as  well  as  deviations  from  the  Imposed  geometrical  constraints.  A  single  sample 
appllcaton  was  presented  In  which  a  known  wing  geometry  was  recovered  for  a  known  target  pressure 
distribution.  The  results  after  six  geometry  Iterations  were  shown.  Convergence  to  the  correct  geometry 
Is  Indicated.  Some  residual  sensitivity  In  the  root  region  of  the  wing  Is  Indicated  In  the  results. 

PAPER  8.  DE  PONTE  presented  an  approach  utilizing  the  ficticious  gas  concept  In  connection  with  a 
field  panel  method  for  compressible  2-0  flows.  The  author  claims  that  the  field  panel  approach  offers  a 
simpler  way  of  Introducing  compressibility  effects.  However,  examples  of  flow  solution  convergence  were 
presented  that  indicated  large  oscillations.  Further,  the  author  Indicated  that  these  oscillations  were 
'FORTRAN  compiler  dependent".  The  author  did  present  an  example  calculation  of  a  2-0  airfoil  modified  to 
be  shock  free  with  the  ficticious  gas  scheme.  While  the  combination  of  ficticious  gas  and  field  panel 
flow  solutions  may  be  unique,  the  work  seemed  rather  sketchy  and  perhaps  reported  prematurely. 


5.0  SESSION  III  -  INVERSE  HETHOOS  /  TURBOMACHINERY,  INTAKES,  DUCTS 


PAPER  3.  VAN  OEN  BRAENBUSSCHE  made  copies  of  his  paper  available  at  the  conference  but  was  unable  to 
be  at  the  conference  for  the  formal  presentation.  A  review  of  the  paper  does  Indicate  a  successful 
approach  for  cascade  blade  design  was  developed.  The  approach  utilizes  a  decoupled  flow 
solver  and  geometry  modifier.  The  geometry  modification  scheme  Is  based  on  a  surface  distribution  of 
vortices  to  drive  an  Initial  velocity  distribution  towards  a  specified  target  distribution.  The  vortex 
distribution  Is  used  to  define  a  normal  flow  which  In  turn  Is  used  In  a  mass  flux  Integral  or  streamline 
slope  Integral  to  define  the  new  geometry  shape.  Examples  are  presented  which  utilize  both  an 
incompressible  potential  flow  solver  and  time-marching  Euler  solver.  Recovery  of  a  known  shape  with  a 
known  velocity  distribution  Is  demonstrated  as  one  of  the  expamples. 

PAPER  10.  CETINKAYA  reviewed  the  development  of  an  Euler  based  oeslgn  method  for  airfoils  and 
cascades.  A  steady,  2-0  Euler  code,  based  on  a  streamline  coordinate  system  and  a  finite  volume 
technique.  Is  used  In  the  method.  The  scheme  Is  seen  to  be  very  similar  to  that  of  Drelal .  Free 
parameters  are  Introduced  Into  the  specified  pressure  distributions  to  account  explicitly  for  LlghthlU's 
constraints.  By  Imposing  a  specified  wall  pressure  distribution  as  a  boundary  condition  In  the  streamline 
calculation,  the  new  surface  geometry  evolves  as  part  of  the  solution  convergence.  Three  examples  are 
presented  which  show  consistent  recovery  of  known  cascade  and  airfoil  geometries  for  known  prescribed 
pressure  distributions.  These  examples  Included  Initial  geometries  exhibiting  embedded  shocks  and  driving 
these  geometries  to  supposedly  shock  free  flows. 
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PAPER  11.  SCHMIDT  described  the  use  of  a  method  based  on  flows  calculated  on  stream  surfaces  of 
revolution.  The  flow  solution  uses  the  full  potential  equation  transformed  Into  the 
potentlal-stream-functlon  plane.  The  discretized  solution  Is  performed  using  successive  line 
overtaxation  and  accelerated  with  a  multigrid  scheme.  Unfortunately,  details  of  the  Inverse  scheme 
Implementation  are  obscure.  Difficulties  relative  to  "111-posedness"  are  mentioned.  The  author  does, 
however.  Indicate  use  of  design  a*  off-design  flow  conditions  as  a  means  of  Improving  off-design 
operational  behaviour.  This  Is  an  approach  to  what  Is  In  reality  a  multi-point  design  problem.  The 
author  also  provides  at  least  a  single  pressure  distribution  correlation  between  a  design  computation  and 
corresponding  experimental  measurement.  This  correlation  does  Indicate  a  good  measure  of  agreement. 

PAPER  12.  JACQUOTTE  presented  an  Inverse  scheme  for  quasi -three-dimensional  flows  through  axial 
cascades.  The  flow  solution  Is  based  on  the  potential  equation  and  Is  solved  using  a  finite  element 
approach.  In  this  scheme  the  Inverse  capability  1$  Implemented  using  the  Dlrlchlet  boundary  condition 
approach  to  Impose  the  specified  pressure  distribution.  The  resulting  solution  normal  flow  to  an  Initial 
surface  Is  Integrated  to  produce  a  surface  displacement  of  the  Initial  geometry.  Both  C-mesh  and  H-mesh 
grid  systems  are  demonstrated.  This  author  also  demonstrates  recovery  of  a  known  geometry.  This 
recovery.  Including  the  convergence.  Is  Illustrated  for  several  different  starting  points.  This  result 
tends  to  Indicate  the  proper  handling  of  Llghthlll's  constraints.  Other  examples  are  also  presented  to 
Indicate  that  apparently  arbitrary  pressure  distributions  did  produce  quite  reasonable  looking  airfoil 
sections.  However,  no  other  results  to  validate  performance  levels  for  the  designed  sections  are 
presented. 

PAPER  13.  BORSES  described  a  three-dimensional  Inverse  method  to  define  a  radial  Inflow  turbine.  The 
method  applies  to  Incompressible  and  Invlscld  flows  and  assumes  that  the  blades  are  Infinitely  thin.  The 
thin  blade  assumption  allows  the  blades  to  be  modelled  by  surface  vortlclty.  The  blade  shape  Is  defined 
by  aligning  the  mean  line  with  the  local  velocity.  Since  this  velocity  depends  on  the  vortlclty 
distribution  the  solution  Is  iterative. 

This  author  was  one  of  the  few  which  utilized  an  Inverse  scheme  to  define  an  aerodynamic  configuration 
which  was  experimentally  validated.  The  author  Is  commended  on  the  thoroughness  of  the  paper  and  the 
completeness  of  the  technical  project.  The  experimental  confirmation  of  an  Improved  Impeller  design  Is  a 
confirmation  of  a  good  computational  approach  to  a  difficult  three-dimensional  problem. 

PAPER  14.  ZANNETTI  presented  a  method  for  the  design  of  three-dimensional  blade  rows  for 
turbomachinery  applications.  The  scheme  utilized  an  Euler  equation  solution  employing  characteristics 
theory.  The  blade  loading  Is  Imposed  In  the  inverse  solution  and  the  blade  camber  Is  calculated. 

6.0  SESSION  IV  -  NUMERICAL  OPTIMIZATION 


PAPER  15.  RIZK  described  an  approach  to  aerodynamic  optimization  In  which  the  Iterative  solution  of 
the  flow  equation  and  the  design  parameter  optimization  are  conducted  simultaneously.  This  rather 
Interesting  approach  avoids  Independent,  and  correspondingly  time-consuming  Iterations  of  the  flow 
solution  and  the  design  parameters.  The  approach  Is  demonstrated  for  a  wind  tunnel  wall  Interference 
study  and  for  a  propeller  design  study.  The  rate  of  convergence  of  the  optimization  problem  appears  to  be 
greatly  enhanced  by  this  simultaneous  approach.  The  examples  presented  Indicate  that  the  optimization  can 
be  accomplished  at  a  cost  of  roughly  11  to  4L  times  the  cost  of  a  regular  analysis  solution  where  L  Is  the 
number  of  design  variables.  The  promise  of  such  an  approach  should  be  pursued  with  further  studies  of  the 
Influence  of  large  numbers  of  design  variables  and  the  effect  of  the  convergence  speed  of  the  underlying 
analysis  method. 

PAPER  16.  VAN  OEN  DAN  presented  a  summary  of  a  spanload  optimization  approach  currently  Implemented 
at  MLR.  The  scheme  represents  a  unified  approach  for  the  preliminary  design  of  multiple  lifting 
surfaces.  Propeller  Influences  on  spanloads  are  also  Included.  Optimization  Is  performed  using  a 
Trefftz-plane  evaluation  of  Induced  drag  plus  form  factor  methods  for  viscous  airfoil  drag.  Minimization 
of  Induced  drag  or  Induced  drag  plus  viscous  drag  can  be  accomplished  subject  to  constraints  on  Items  such 
as  lift,  moment,  and  airfoil  characteristics.  The  method  Is  demonstrated  for  numerous  configurations 
Including  canard  and  three-surface  configurations,  and  wing  tall  variations.  The  method  Is  clearly  a 
valuable  preliminary  design  tool.  It  Is  an  accumulation  of  a  number  of  classical  procedures  that  many 
aircraft  companies  utilize  in  one  manner  or  another.  The  present  approach  would  appear  to  be  the  result 
of  a  conscious  effort  to  cast  these  preliminary  design  procedures  In  a  consistent  optimization  format  with 
considerable  Flexibility. 

PAPER  17.  VAN  EGMOND  described  the  development  and  application  of  parameterized  airfoil  pressure 
distributions  that  serve  as  target  pressure  distributions  for  Inverse  design  studies.  This  approach  has 
been  developed  to  help  systematically  answer  the  question  as  to  what  Is  the  correct  pressure  distribution 
to  specify  In  an  Inverse  design  problem.  The  parametric  representation  of  an  airfoil  pressure 
distribution  Is  reviewed.  Several  means  of  numerical  optimization  are  provided.  The  method  consists  of 
the  parametric  pressure  distribution  representation,  boundary  layer  calculations  for  viscous  drag,  a 
simple  formula  for  approximate  shock  drag,  and  an  approximate  means  for  constraining  airfoil  thickness 
based  on  a  pressure  distribution  Integral.  Since  these  Ingredients  are  relatively  simple  and  Inexpensive 
In  terms  of  computational  resources,  many  Iteratons  can  easily  be  accomodated  In  the  optimization  scheme. 

The  typical  use  of  this  method  would  seem  to  be  initiated  with  an  analysis  solution  on  an  existing 
airfoil.  A  parametric  best  fit  of  the  resulting  pressure  distribution  Is  computed.  This  best  fit 
distribution  Is  then  numerically  Improved  ("optimized”)  using  the  elements  described  above.  The  final 
pressure  distribution  Is  then  used  as  an  Input  target  pressure  distribution  using  an  Inverse  design 
method.  Example  applications  to  laminar  flow  airfoils,  a  llebeck  high-lift  airfoil,  and  a  transonic 
airfoil  are  described.  Only  calculations  are  presented;  no  experimental  verification  Is  demonstrated.  In 
fact,  the  differences  In  lift  and  drag  shown  for  the  optimized  airfoil  configurations  are  small  and  may  be 
within  the  computational  and  experimental  error  band. 


This  scheme  begins  to  resemble  an  approach  In  which  an  Inverse  airfoil  code  is  Iteratively  driven  by  a 
parametric  pressure  distribution  which  Is  Itself  driven  simultaneously  by  a  numerical  optimization 
scheme.  Such  an  approach  seems  highly  desirable  to  the  reviewer.  However,  the  MLR  scheme  presented  by 
Van  Egmond  does  MOT  drive  an  Inverse  code  during  the  numerical  otlmlzatlon.  Instead  It  makes  use  of  the 
approximate  elements  described  above  and  upon  completion  of  the  optimization  makes  use  of  the  Inverse 
solution  only  once.  While  the  MLR  approach  Is  clearly  an  economical  one,  the  additional  aproxlmatlons 
Introduced  may  limit  Its  utility  In  cases  of  small  refinements. 

PAPER  18.  6HIELMI  provided  a  description  of  a  numerical  optimization  scheme  for  airfoils  at  multiple 
operating  points.  The  scheme  utilized  the  well  known  COMMIN  algorithms  as  the  design  parameter  driver. 
Explicit  algebraic  shape  functions  are  used  to  perturb  the  airfoil  contours.  Specific  reference  Is  made 
to  the  emphasis  on  the  use  of  constraints  rather  than  a  complicated  objective  function.  The  example 

presented  was  an  airfoil  design  problem  for  a  military  trainer  with  two  transonic  design  points.  The  flow 
solver  used  In  this  effort  was  the  two-dimensional  potential  solver,  FL06,  of  Jameson.  Additionally,  the 
author  in  the  written  paper  makes  mention  of  the  use  of  expert  system  programming  to  effectively  implement 
the  scheme  and  retain  man-ln-loop  capability. 

PAPER  19.  RENEAUX  described  an  approach  similar  to  the  previous  paper.  The  ONERA  scheme  Is  based  on 
the  application  of  CONMIN  combined  with  the  two-dimensional  airfoil  analysis  method  of  Garabedlan  and 
Korn.  The  numerical  optimization  scheme  was  applied  to  a  helicopter  rotor  design  problem.  The  scheme  was 
first  used  to  develop  Improved  airfoil  sections  using  two  design  points.  One  point  corresponds  to  an 
advancing  blade  condition,  while  the  second  point  corresponds  to  a  retreating  blade  condition.  The  blade 
section  Improvements  that  were  computationally  developed  were  experimentally  verified  by  wind  tunnel 
test.  The  Improvements  are  significant  enough  to  be  well  beyond  the  computational  and  experimental  error 
band.  Consequently,  the  Investigators  had  a  high  probability  of  success.  The  authors  are  to  be  commended 
for  the  completeness  of  the  airfoil  optimization  effort. 

A  second  numerical  optimization  application  Is  also  demonstrated.  The  second  application  uses  the 
same  optimizer  coupled  to  a  rotor  performance  method  based  on  blade  element  theory.  It  Is  emphasized  that 
such  an  approach  Is  attractive  due  to  the  many  flight  conditions  and  constraints  which  must  be  observed  In 
the  rotor  design  process. 

PAPER  20.  BOCK  presented  a  summary  of  numerical  optimization  cases  Including  a  supersonic  airfoil 
section,  a  supersonic  body  of  revolution,  a  transonic  airfoil,  and  a  subsonic  multi-element  high-lift 
airfoil.  The  paper  appears  to  be  a  thorough  documentation  of  these  examples  and  Indicates  attention  being 
paid  to  real  design  problems.  Clearly,  the  author  has  found  that  a  numerical  optimization  approach  offers 
value  to  the  aerodynamic  designer. 

The  author  did,  however,  Indicate  several  points  which  deserve  some  attention.  First,  an  Initial 
comparison  Is  made  between  gradient  methods  and  evolution  methods  for  numerical  optimizers.  The 
conclusion  Is  reached  that  the  gradient  method  Is  far  superior  as  a  result  of  the  extended  number  of 
Iterations  required  by  the  evolution  method.  While  the  conclusion  may  be  correct  for  specific  problems, 
this  Implementaton  of  the  evolution  method  would  appear  to  be  such  that  an  Inefficient  search  Is 
accomplished.  The  evolution  method  should  be  a  trivial  scheme  to  Implement  and  has  been  shown  to  converge 
rapidly  In  other  Investigations.  Second,  In  the  transonic  airfoil  case  a  numerical  optimization  scheme  Is 
presented  In  which  the  design  variables  Include  an  upper  surface  bump  function  defined  by  a  ficticious  gas 
scheme  to  drive  towards  a  shock-free  airfoil  contour.  Such  an  approach  Is  believed  by  the  reviewer  to  be 
overly  constraining.  The  ficticious  gas  scheme  Is  an  Interesting  academic  exercise  but  does  not  provide 
design  capability  for  the  entire  airfoil  surface.  Instead,  It  focusses  on  the  supersonic  bubble  and  the 
attainment  of  shock-free  flow.  Airfoils  designed  for  shock-free  flows  are  notorious  for  poor  off-design 
behaviour.  Consequently,  a  better  optimization  can  be  accomplished  If  the  scheme  addresses  the  Issue  of 
shocked  flows  at  multiple  design  points.  Third,  an  observation  Is  made  that  convergence  of  the  objective 
function  Is  not  clear  for  several  of  the  cases  presented. 

PAPER  21.  OESTARAC  reviewed  examples  of  numerical  optimization  for  transonic  wing  design  problems. 
Preliminary  efforts  relating  two-dimensional  and  three-dimensional  results  provide  an  Indication  of  the 
background  utilized  In  the  wing  examples  presented.  The  three-dimensional  wing  optimizations  presented 
Include  the  development  of  the  Inboard  area  of  a  three-dimensional  swept  wing  and  the  development  of  local 
treatments  required  to  Integrate  wing  mounted  nacelle/pylon  configurations.  These  two  problems  are 
classical  wing  aerodynamic  design  problems.  The  results  presented  for  numerically  optimized  configuration 
variations  Indicate  a  rapid  means  to  develop  Improved  geometry  solutions  for  these  problems.  The 
convergence  results  presented  In  this  paper  Indicate  significant  convergence  Is  Indeed  being  obtained  by 
the  numerical  optimization  scheme. 

PAPER  23.  HUDDLESTON  presented  the  progress  on  a  rather  unique  application  of  numerical  optimization 
coupled  to  a  CFD  method.  This  appllcaton  Is  targeted  at  the  development  of  a  simulation  for  engine  Inlet 
testing  In  a  ground-based  facility.  The  simulation  Is  to  provide  an  Inlet  onset  flowfleld  that  best 
matches  that  which  a  fighter  vehicle  forebody  can  create.  One  of  the  unique  aspects  of  this  work  Is  that 
the  optimization  objective  function  Is  based  on  off-body  flow  variables.  A  least-squares  type  of 
evaluation  Is  used  to  make  the  evaluation  of  the  objective  function  consistent  with  the  design  variables 
Implemented.  Three  test  cases  are  reviewed  In  which  the  off-body,  least-squares  approach  Is  used  to 
recover  a  known  and  consistently  generated  target  solution.  The  scheme  demonstrates  rapid  convergence  for 
these  cases.  However,  the  author  does  Indicate  some  pre-conditioning  had  to  be  utilized  In  cases 
presented.  The  application  of  this  approach  to  the  real  design  problem  will  be  quite  Interesting. 
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7.0  GENERAL  COMMENTS 


This  Technical  Specialists'  meeting  on  aerodynamic  design  and  optimization  methods  was  timely  and  well 
developed.  The  program  Included  papers  which  dealt  with  most  approaches  to  aerodynamic  design  and 
optimization.  The  subject  presentations  Included  what  could  be  termed  classical  approaches  to  Inverse 
aerodynamic  methods  as  well  as  newly  developed  and  highly  promising  schemes.  The  program  contributions 
were  also  from  a  broad  International  spectrum  and  covered  a  wide  variety  of  applications  from  the 
classical  airfoil  problem  to  complex  turbomachinery  problems.  This  variety  Indicates  the  broad  base  of 
such  applications  currently  being  exercised  throughout  the  global  aerospace  community. 

Clearly,  the  quality  of  the  studies  varied  and  the  quality  of  the  design  results  also  varied.  Some 
authors  only  reported  preliminary  results  from  analytical  efforts.  If  the  design  or  optimization  results 
are  modest  changes  from  a  base  conf Iguraton,  It  Is  often  hard  to  prove  that  the  computational 
conf Iguratlon  refinement  Is  more  than  bouncing  around  In  the  computational  and  experimental  error  band. 
Others  reported  thorough  analytical  results  backed  up  by  experimental  confirmation  of  the  design 
aerodynamics.  This  latter  completeness  Is  to  be  commended. 

It  Is  obvious  that  computational  design  and  optimization  efforts  will  continue  to  grow  as  pressures 
for  Increased  automation  and  design  capability  are  felt.  Increased  global  competition  will  continue  to 
add  such  pressures  and  the  design  aerodynamlclst  community  will  Increasingly  rely  on  such  methods  to 
Improve  design  efficiency.  It  1<  recommended  that  AGARO  continue  Its  Interest  In  such  methods  and 
formulate  future  events  for  timely  dissemination  of  results  on  an  International  scale. 
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